Identification of Dark Matter particles with LHC and direct detection data 
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Dark matter (DM) is currently searched for with a variety of detection strategies. Accelerator 
searches are particularly promising, but even if Weakly Interacting Massive Particles (WIMPs) are 
found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the 
bulk of the DM in the Universe $1dm- We show that a significantly better reconstruction of the DM 
properties can be obtained with a combined analysis of LHC and direct detection (DD) data, by 
making a simple Ansatz on the WIMP local density p^o , i.e., by assuming that the local densiy scales 
with the cosmological relic abundance, (pjo/pdm) = (fi^o /Hdm)- We demonstrate this method in 
an explicit example in the context of a 24-parameter supersymmetric model, with a neutralino LSP 
in the stau co-annihilation region. Our results show that future ton-scale DD experiments will allow 
to break degeneracies in the SUSY parameter space and achieve a significantly better reconstruction 
of the neutralino composition and its relic density than with LHC data alone. 
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I. INTRODUCTION 

Identifying the nature of the dark matter (DM) re- 
mains one of the central unsolved problems in modern 
particle physics and cosmology. A generic Weakly In- 
teracting Massive Particle (WIMP) is among the best- 
motivated possibilities since it can be thermally pro- 
duced in the early Universe in the right amount to ac- 
count for the observed DM density. Indeed, many theo- 
ries for Physics beyond the Standard Model contain vi- 
able WIMP candidates, as is the case of Supersymmetry 
(SUSY) when the lightest SUSY particle (LSP) is the 
lightest neutralino (a linear superposition of the super- 
symmetric partners of the gauge and Higgs bosons) [l|-[l] . 

DM can be searched for in various ways. One possi- 
bility is attempting a direct detection, through its scat- 
tering off nuclei inside an underground detector. Many 
experiments have been running or arc under construc- 
tion which are mostly sensitive to the spin-independent 



part of the WIMP-nucleus cross section, a? 



Among 



these, the DAMA/LIBRA collaboration reported a pos- 
sible DM signal [a, Q. However, its interpretation in 
terms of the elastic scattering of a WIMP with a mass 
around 10 - 100 GeV and a^_ p ~ 10~ 3 - 10~ 5 pb has 
been challenged by other experiments, such as the Co- 
GeNT @,1, CDMS @, 03, XENON [n[ and ZEPLIN 
fl2| . The CoGcNT collaboration has itself recently re- 
ported an irreducible excess of low-energy events which 
could also be understood as due to the scattering of a very 
light WIMP (see also Ref. [HJ]), but this intepretation 
has in turn be put under pressure by the XENON-100 
results, obtained with a fiducial target mass of 40 kg and 
1 1 days of exposure [14[ . Finally, the recent results from 
the CDMS-II collaboration show two events compatible 
with a WIMP signal, although these results are still sta- 



tistically inconclusive [10| . 

The future increase of the sensitivity may clarify the 
situation, but it is becoming clear that several indepen- 
dent pieces of evidence will be necessary to claim discov- 
ery of DM. In fact, even if in principle the WIMP mass 
and scattering cross section can be determined with some 
accuracy after its direct detection in one direct detec- 
tion experiment, provided that the measured event rate 
is large and the WIMP mass is small [15|, [Ty] , a second 
direct detection with a different target would actually 
allow a much more precise determination of the WIMP 



17 1 , and if the new target is sensitive to the spin- 



dependent contribution of the WIMP-nucleus cross sec- 
tion it could even be used to discriminate among WIMP 
candidates [LSI ]. 

Another possibility consists in looking for the products 
of DM annihilation (e.g., high energy neutrinos, gamma- 
rays or antimatter) and thus indirectly reveal the presence 
of the DM 0, [j] . We leave the discussion of this search 
strategy to a forthcoming work, where we will present the 
constraints that can be set on the DM parameter space 
from the observation (or non-observation, see also fl9l ]) 
of DM annihilation radiation [2fJ . 

Finally, collider experiments, most notably the Large 
Hadron Collider (LHC), will explore the nature of 
Physics at the TeV scale, where many of the extensions 
of the SM that propose DM candidates would manifest 
themselves. The detection of new Physics in particle col- 
liders can provide crucial information about DM. For ex- 
ample, the mass and spin of the LSP could be deter- 
mined through the study of kinematic variables |2l], [22| ■ 
However, to prove that the newly discovered particles ac- 
count for all (or most) of the DM in the Universe, is a 
challenging task. In fact, although particle accelerators 
can provide some information about the neutralino relic 



density [23[, it was found that in many cases the LHC 
would be unable to determine the precise composition of 
the neutralino, leading to an unreliable prediction of its 
relic abundance or to the occurrence of multiple solutions 
spanning several orders of magnitude, thus not allowing 
to establish whether or not it is the DM (see also Ref. 
[24j and references therein). 

One possibility is to build a new collider, such as the 
proposed International Linear Collider (ILC), that would 
allow a much more precise evaluation of the supersym- 
mctric masses and couplings, and a better determination 
of the inferred relic density, as argued by the authors of 
Ref. |23[ . However, this machine will not be available in 
the near future, and it is therefore crucial to devise strate- 
gies that can be implemented as soon as new particles 
are discovered at the LHC. Fortunately, direct detection 
experiments are expected to greatly improve their sensi- 
tivity in the next few years and start probing interesting 
regions of the supcrsymmetric parameter space. In case 
of discovery, it will certainly be reassuring if the mass re- 
constructed from direct detection experiments matched 
the value obtained from accelerator measurements, since 
it would prove the existence of a particle which is sta- 
ble over cosmological timcscalcs. The error on the mass 
reconstructed from direct detection experiments depends 
on the DM particle parameters, and on the experimental 
setup, and the interested reader can find a detailed anal- 
ysis in Refs.fla. Ha |. But one can do much more than 
checking the compatibility of the two mass determina- 
tions. We show here that a combined analysis of the two 
data sets will allow a much better reconstruction of the 
DM properties, and a convincing identification of DM 
particles. 

Although the strategy discussed here is model- 
independent, we work out an explicit example in the 
context of a 24-parameters supcrsymmetric model, with 
a neutralino LSP in the stau co-annihilation region. 



II. THEORETICAL FRAMEWORK AND LHC 
DATA 

We work within the framework of the minimal super- 
symmetric extension of the Standard Model (MSSM), for 
which we adopt a low energy parametrization in terms of 
24 parameters, corresponding to its CP-conserving ver- 
sion. The input parameters are the coefficients of the tri- 
linear terms for the three generations, the mass terms for 
gauginos (for which no universality assumption is made) , 
right-handed and left-handed squarks and leptons, the 
mass of the pseudoscalar Higgs, the Higgsino mass pa- 
rameter /i, and finally the ratio between the vacuum ex- 
pectation values of the two Higgs bosons tan/3. 

If searches for new Physics at the LHC are consistent 
with a SUSY scenario, the study of different kinematical 
variables will allow us to determine some properties of 
the SUSY spectrum. In particular, the masses of several 
particles or mass-splittings between them could be ex- 



tracted, with a precision that obviously depends on the 
properties of the specific point of the parameter space. 
These measurements can then be used as constraints on 
the 24-dimensional SUSY model, in order to determine 
the regions of the MSSM parameter space which are con- 
sistent with such a measurement. This can be done by 
applying B ayes' theorem 



p(x|d) 



p(d|x)p(x) 
p(d) 



(1) 



which updates the so-called prior probability density 
p(x), encapsulating the knowledge of the 24-dimcnsional 
space before taking into account the experimental con- 
straints, d, into the posterior probability function (pdf) 
p(x|d). The latter describes the probability density as- 
signed to a generic 24-dimensional point x once the data 
have been taken into account via the likelihood function 
p(d|x). Furthermore, on the RHS of Eq. (JTJ, p(d) is the 
Bayesian evidence which, in our case, can be dropped 
since it simply plays the role of a normalization constant 
for the posterior in this context (see [25[ for further de- 
tails). 

The marginal pdf of a particular subset (as e.g. only 
one) of the 24 parameters defining x can be obtained by 
integrating over the remaining directions: 



p(x*|d) = f p(x\d)dx 1 ...dx i - 1 dx i+1 . 



.dx 2 



(2) 



The posterior encodes both the information contained 
in the priors and in the experimental constraints, but, 
ideally, it should be largely independent of the choice of 
priors, so that the posterior inference is dominated by 
the data contained in the likelihood. If some residual de- 
pendence on the prior p(x) remains this should be con- 
sidered as a sign that the experimental data employed 
are not constraining enough to override completely dif- 
ferent plausible prior choices and therefore the result- 
ing posterior should be interpreted with some care, as it 
might depend on the prior assumptions. The probabil- 
ity distribution for any observable that is a function of 
the 24 SUSY parameters /(x) can also be obtained since 
p(/|d) = <5(/-/(x))p(x|d). 

For the practical implementation of the Bayesian anal- 
ysis sketched above we employed the SuperBayeS code 
[26(, extending the publicly available version 1.35 to han- 
dle the 24 dimensions of our SUSY parameter space. To 
scan in an efficient way the SUS Y p arameter space we 
have upgraded the MultiNest [23, [28[ algorithm included 
in SuperBayeS to the latest MultiNest release (v 2.7). 
MultiNest is a multi-modal implementation of the nested 
sampling algorithm, which is used to produce a list of 
samples in parameter space whose density is proportional 
to the posterior pdf of Eq. ([T]). For further information 
on nested sampling we refer the reader to the appendix 
of Ref. [29[ and references therein. 

For the present work we have chosen a specific bench- 
mark point in the MSSM parameter space, corresponding 



to the low-energy extrapolation of model LCC3 defined 
in Ref. [23J]. This benchmark is representative of SUSY 
models in the co-annihilation region, where the lightest 
neutralino is almost degenerate in mass with the light- 
est stau. In this region, co-annihilation effects reduce 
the neutralino relic abundance down to values compati- 
ble with the results from the WMAP satellite [H, and 
therefore, the mass difference between the neutralino and 
the lightest stau is a fundamental parameter for the re- 
construction of the relic density. It has been shown [23[ 
that for this benchmark point LHC would be able to pro- 
vide a measurement of the masses of a good part of the 
SUSY spectrum, including the two lightest neutralinos 
(see Ref. |31[ for an extension of this analysis to the case 
of the ILC). However the masses of some particles (most 
notably the two heaviest neutralinos and both charginos) 
would not be measured. The set of measurements that 
we use as constraints in our analysis corresponds to that 
in Table 6 of Ref. [23j l , which assumes an integrated 
luminosity of 300 fb -1 . Furthermore, as pointed out in 
Ref. [321 ] , the neutralino-stau mass difference can be mea- 
sured with an accuracy of 20% with 10 fb -1 luminos- 
ity in models where the squark masses are much larger 
than those of the lightest chargino and second-lightest 
neutralino, as is our case. We therefore also include a 
measurement of the neutralino-stau mass difference in 
our likelihood. For convenience, we summarize in Table 
U the set of LHC measurements on which we build our 
likelihood. Each of the constraints listed in Table Q] is 
implemented in the likelihood as an independent Gaus- 
sian distributed measurement around the true value p, for 
that observable, with standard deviation er, as given in 
Table H 

Compared with previous Bayesian studies in which 
only precision tests of the Standard Model are considered 
as experimental constraints [23, I33I435J , we are assuming 
here a scenario in which LHC reports a quite stringent 
collection of measurements. For this reason our poste- 
rior constraints are quite tight and we expect the prior 
dependence of our results to be very mild. This is con- 
firmed by the inspection of the profile likelihood, which 
agrees well with the posterior pdf (see [29| for a detailed 
discussion). This indicates that volume effects from the 
prior are unlikely to be playing a major role given the 
strong constraints we assume for our benchmark point. 



III. FUTURE DIRECT DETECTION DATA 

In the simulation of a direct detection experiment we 
assume a future signal giving a WIMP detection, namely 



Mass 


Benchmark value, /j, 


LHC error, a 


m(Xi) 


139.3 


14.0 


m(xS) 


269.4 


41.0 


m(ei) 


257.3 


50.0 


m(pi) 


257.2 


50.0 


m(h) 


118.50 


0.25 


m{A) 


432.4 


1.5 


m(ri) - 


m(x?) 16.4 


2.0 


m(u R ) 


859.4 


78.0 


m(d R ) 


882.5 


78.0 


m(s R ) 


882.5 


78.0 


m(c R ) 


859.4 


78.0 


m(uL) 


876.6 


121.0 


m(d,L) 


884.6 


121.0 


m(J L ) 


884.6 


121.0 


m(c L ) 


876.6 


121.0 


m(bi) 


745.1 


35.0 


m(b 2 ) 


800.7 


74.0 


m(ti) 


624.9 


315.0 


m(g) 


894.6 


171.0 


m(e 2 ) 


328.9 


50.0 


mip^) 


328.8 


50.0 



TABLE I: Sparticle spectrum (in GeV) for our benchmark 
SUSY point and relative estimated measurements errors at 
the LHC (standard deviation er). 



1 The exact values that we are using arc slightly different from 
those in Table 6 of Ref. |2.'1 since we are not deriving the mass 
spectrum from the low energy extrapolation of a cMSSM point. 
On the contrary our reference SUSY model is defined at low 
energy to be near LCC3. 



a certain number of events N and a corresponding set 
of recoil energies {-Ei}i=i,...,jv- The total number N of 
simulated events is the sum of both background events 
(mainly interactions of detector nuclei with neutrons 
from surrounding rock, from residual contaminants or 
from spallation of cosmic muons) and recoils due to DM. 
For concretcness, we will exemplify the method in the 
case of an experiment akin to the 1-ton scale Supcr- 
CDMS experiment j36[. We simulated the differential 
number of background events as in Ref. [37|. Since the 
capability of a simulated direct detection experiment to 
reconstruct the DM properties (see Refs. yJI, [l|| |37( for 
more details) is worse in the case of a constant back- 
ground distribution than for an exponential one, we only 
consider the case of energy-independent background re- 
coil spectrum in order to be conservative. Therefore, 
we adopt a constant background differential spectrum 
(dNback/dE) = const which is normalized so that, when 
binning the spectrum in 9 bins of 10 kcV width (from 
JSth = 10 keV to -Emax = 100 keV) the number of back- 
ground events in the first bin is the same as the number 
of DM signal events there. 

The expected number of events A for our benchmark 
model and for an exposure e = 300 ton day is obtained 
by integrating the sum of the differential rate of WIMP 



Target 


A 


e 


£th 


i'max 


Px 


A 


Ge 


73 


300 ton day 


lOkeV 


100 keV 


0.385 GeV cm" 3 


638 



TABLE II: Relevant quantities for a SuperCDMS-like direct 
detection experiment. The quantity A gives the expected 
number of WIMP recoils for our SUSY benchmark model. 



and background events 



-k'th 



dR v dR 



hack 



dE 



dE 



dE. 



(3) 



The dependency of the WIMP event rate on the phys- 
ical quantities in the problem becomes apparent in the 
following parametrization |38f 



where 



and 



dR x 
~dE 



Ro 



c 1 Roe- E ^ Eoc ^F 2 (E), 



a s x 1 _ p p x A 2 c 2 (m x +m p f 
y/%m 3 m,pVo 



En 



2m 2 ,v 2 Am p 
(m x + Am p ) 2 c 2 



(4) 



(5) 



(6) 



Here, p x is the local WIMP density, A is the mass number 
of the target nuclei (A = 73 in the case of Germanium), 
m p is the proton mass, vo is the characteristic WIMP 
velocity and F 2 (E) denotes the nuclear form factor. A 
discussion on the values of the parameters C\ and c-i and 
the functional form of F(E) can be found in Refs. [15|, 
Ho . l3q . The specific values of the quantities for our case 
study are summarized in Table [Til 

In order to combine the result of a direct detection ex- 
periment with LHC data, we run an additional scan of the 
SUSY parameter space including in the likelihood func- 
tion an additional Poisson-distributed term that com- 
pares the number of events and their spectral shape pre- 
dicted in each point in parameter space with the recoil 
spectrum corresponding to the benchmark value of Table 
HIl The overall background rate and its spectral shape are 
assumed to be known. 

As shown by Eqs. ||3J|-||SJ|, the number of detected 
events is proportional to the product of the WIMP- 
proton cross section and the local DM density A ex 
&x-pPx' Therefore, unless one specifies the value of p x , 
any information on the number of events leaves the scat- 
tering cross section practically unconstrained. 

We propose two different strategics to specify p x : 



1. Consistency check: we impose that 



Px — Pdm , 



(7) 



and we adopt for this quantity the value obtained 
in a recent paper by Catena and Ullio [39| , through 



a careful analysis of dynamical observables in the 
Galaxy, namely p x = 0.385 GeV cm -3 (see also 
[40l - [43l | ) . Although this assumption completely re- 
moves the degeneracy between cr?} and p x , it 
forces the identification of neutralino with the DM 
particle, irrespectively of the value of its thermal 
relic density. This is therefore equivalent to as- 
suming that, a non-standard cosmological history 
of the Universe can correct any excess or deficit in 
the thermal relic density and make it agree with 
the WMAP result, for example, either by invoking 
late injection of entropy, non-thermal production 
through late-decaying particles (such as a modulus 
or a gravitino |44l0. scenarios with a low- reheating 
temperature Jia (see also Ref. [46] ) or a faster ex- 
pansion rate [47|, |48| . For these reasons this Ansatz 
must be considered as a consistency check rather 
than a proof of the identification of DM particles. 

2. Scaling Ansatz: we assume that the local density of 
the neutralino scales with the cosmological abun- 
dance. More precisely, we propose the following 
Ansatz 



p x o/p 



DM 



= a 



„/n 



DM- 



(8) 



This Ansatz is strictly valid in the reasonable case 
where the distribution of neutralinos in large struc- 
tures, and in particular in the Galaxy, traces the 
cosmological distribution of the DM. This Ansatz 
is obviously true if neutralinos contribute all the 
DM in the Universe, but is also valid in the case 
where the neutralino is a subdominant component 
of DM, provided that DM behaves, as expected, as 
a cold collisionlcss particle. As shown below, this 
simple assumption is powerful tool to remove de- 
generacies in the parameter space. 

The reconstruction of the neutralino relic density is 
shown in Fig.Q] The left panel corresponds to the case 
where only LHC constraints are considered. Consistently 
with previous analyses [231 ] . multiple peaks can be ob- 
served, as a consequence of degeneracies in the SUSY 
parameters space that the LHC constraints are unable 
to break. In particular, the two observed peaks corre- 
spond to neutralinos with different composition: mostly 
Wino and mostly Bino, from left to right. This is a con- 
sequence of the fact that the LHC is assumed to be able 
to measure only the two lightest neutralino states, but 
not the two more massive ones or the charginos. The 
true value of the relic density for our benchmark point 
{Vt x oh 2 = 0.176), represented by a diamond in Fig. [TJ is 
indeed inside the peak corresponding to mostly Bino dark 
matter. Although this value is about 60% larger than the 
relic abundance measured by the WMAP satellite j30j |. 
we expect our results to remain qualitatively correct for 
other points in the co-annihilation region leading to the 
correct cosmological relic abundance. As commented 
above and already pointed out in previous works [23j the 
better reconstruction of sparticle masses at the ILC could 
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FIG. 1: Pdf of the neutralino relic density obtained from LHC data only (left panel), LHC plus direct detection data for a fixed 
local density (central panel) and LHC plus direct detection data under the scaling Ansatz. The best fit point is shown by the 
encircled black cross, while the true value is given by the yellow/red diamond. 
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FIG. 2: As in Fig. Q] but in the at 
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plane. From left to right: constraints using LHC data alone, adding direct 



detection data assuming a fixed local DM density, and with our scaling Ansatz for the local density. The inner and outer 
contours enclose 68% and 95% probability regions, respectively. The probability distributions have been smoothed with a 
Gaussian kernel for display purposes. The best fit point is shown by the encircled black cross, while the true value is given by 
the yellow/red diamond. In the center and right panels we show for reference the LHC-only contours of the left panel (light 
grey), along with the directions along which the different Ansatze break the degeneracy (dashed lines). 



allow a more precise determination of the neutralino relic 
abundance, potentially removing some of these degenera- 
cies. However, this information would only be available 
after a much longer period of time. 

The constraints from LHC only data are also shown 
in the left panel of Fig. [2j in the plane c?, 1 vs il^oh 2 , 
where the true value of those quantities is given by 
{il^oh 2 = 0.176, ct| : _ p = 7.1 x 1CT 8 pb). The left- 
most region corresponds to a neutralino which has a lead- 
ing Wino component, thereby displaying a smaller relic 
abundance, whereas the region towards larger relic abun- 
dance corresponds to Bino-likc ncutralinos, for which the 
scattering cross section is also slightly smaller. 

In the central and right panels of the two figures, we 
show the impact of adding information from direct de- 
tection experiments. These plots have been obtained by 
statistical posterior re-sampling of the LHC only scan, 



adding the relevant Ansatze and the likelihood function 
of a direct detection experiment as specified above. The 
central panels correspond to the assumption 1, or consis- 
tency check. This amounts to fixing the local neutralino 
density, and therefore we expect that only regions along 
a direction of constant c?} to survive after direct de- 
tection data are implemented. This can be understood 
as follows: for a given number of measured events, and 
a fixed local density, there is only a range of values of 
ctS : _ p that are compatible with the measurement. Notice 
that, as explained above, this Ansatz does not further 
constrain the neutralino thermal relic density. In this 
case the pdf for the neutralino relic density still displays 
the two maxima, corresponding to the two peaks in Fig. [I] 
and the two "islands" in Fig.[5J This is due to the fact 
that the neutralino can have a similar scattering cross 
section for both compositions and therefore (if the fact 



that it might be a subdominant DM component is not 
properly taken into account) could account for the same 
detected rate. 

The most interesting case is the one that corresponds 
to our assumption 2, namely the scaling Ansatz, which 
represents the most important result of this paper. When 
the appropriate scaling of the local density is applied, the 
Ansatz cuts the parameter space along a direction 



°£,«n«. 



(9) 



due to the fact that for a fixed number of events <rS oc 



p x l and that under the scaling Ansatz p x oc 



f2-po. 



The 



dramatic consequences of this simple Ansatz are shown 
in the right-most panels of both figures. Models corre- 
sponding to a low relic density are essentially ruled out, 
because under the scaling Ansatz they correspond to a 
low local density. Given a number of observed events 
in direct detection searches, a low local density would 
require a larger scattering cross section, which is incom- 
patible with LHC constraints. As a consequence, the pa- 
rameter space region corresponding to a neutralino that 
is mostly Wino can now be ruled out with high confi- 
dence, thereby leading to a much better reconstruction 
of the DM composition than it would be possible under 
the consistency check Ansatz. 

We note that if the reconstructed relic density matches 
the observational determination of Hdm, this procedure 
also validates the standard cosmological history, and con- 
strains deviations from the standard expansion rate at 
the epoch of DM freeze-out. Conversely, a mismatch 
between the reconstructed relic density and Hdm would 
point towards a multi-component DM sector, or a non- 
standard expansion rate (see e.g. Ref. [49J and refer- 
ences therein). 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated the effect of combining informa- 
tion from accelerator searches with data from direct DM 
detection, assuming realistic measurements at the LHC 
and in a Germanium detector with an exposure of 300 
ton days. 

An interesting question is whether the systematic and 
statistical errors on this quantity and on other relevant 
physical quantities entering in the calculation of the event 
rate for direct detection experiments can spoil the re- 
construction procedure presented here. For instance, we 
have assumed a Maxwcllian distribution for the velocity 
dispersion of DM particles, but a more refined analysis 
should keep into account the uncertainties on this quan- 
tity. Fortunately, recent estimates based on numerical 
simulations, suggest that the small measured deviations 
from a Maxwcll-Boltzmann distribution lead to errors of 



10% or less on the recoil rate, and they are therefore sub- 
dominant with respect to other uncertainties, such as the 
error on the nuclear form factor [5(j], and especially the 
error on the observed DM local density. The most impor- 
tant effect of such uncertainties, once marginalized over, 
would be to widen the pdf's of Fig.[2]in the vertical direc- 
tion by less than 10%, if one considers only the statistical 
error on /3dm derived in Ref. [33 ] , and by up to a factor 
of two if one also considers the systematic error due to 
halo triaxiality [43| ■ Since the vertical thickness of the 
contours in the 2D posterior of Fig. [2] is approximately 
equal to a factor of 2, we expect that including these un- 
certainties would not modify qualitatively the marginal 
posterior distribution for H, x oh 2 , and our results would 
still apply. A more detailed discussion of these effects 
is beyond the scope of this paper, and we leave it to a 
separate upcoming work. However we explicitly studied 
the effect of varying the value for the mass of the top, 
including it as a nuisance parameter in the likelihood. 
The variation in the reconstructed pdf for the neutralino 
relic abundance and neutralino-nucleon scattering cross 
section is negligible. 

We stress once more the importance of combining dif- 
ferent types of experiments. The specific case discussed 
here shows that when reasonable assumptions are made 
to link the local density to the relic abundance, a com- 
bined analysis of data from accelerators and direct detec- 
tion experiments allows a significantly better reconstruc- 
tion of the DM properties. 

This is true in the co-annihilation region discussed 
here, but it will provide important information for any 
SUSY scenario, and more in general for any new physics 
scenario. Even in cases where the LHC data are sufficient 
to pinpoint the underlying DM scenario, direct detection 
experiments can corroborate the results, and they can 
also be used to identify deviations from the standard ex- 
pansion rate of the Universe at freeze-out that would ap- 
pear as an inconsistency between the f2^o inferred from 
LHC data and cosmological measurements. 
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